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Abstract 

We introduce a new family of indecomposable positive linear maps based on 
entangled quantum states. Central to our construction is the notion of an 
unextendible product basis. The construction lets us create indecomposable 
positive linear maps in matrix algebras of arbitrary high dimension. 
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I. INTRODUCTION 



One of the central problems in the emergent field of quantum information theory [jTJ is 
the classification and characterization of the entanglement (to be defined in section |TT|) of 
quantum states. Entangled quantum states have been shown to be valuable resources in 
(quantum) communication and computation protocols. In this context it has been shown 
that there exists a strong connection between the classification of the entanglement of 
quantum states and the structure of positive linear maps. Very little is known about the 
structure of positive linear maps even on low dimensional matrix algebras, in particular the 
structure of indecomposable positive linear maps. We denote the n x n matrix algebra as 
M„(C). The first example of an indecomposable positive linear map in M^(C) was found 
by Choi ||. There have been only several other examples of indecomposable positive linear 
maps (see Ref. |3J for some recent literature); they seem to be hard to find and no general 
construction method is available. In this paper we make use of the connection with quantum 
states to develop a method to create indecomposable positive linear maps on matrix algebras 
M„(C) for any n > 2. Central in this construction is the notion of an unextendible product 
basis, of which there exist examples in arbitrary high dimensions |5|J|. Unextendible product 
bases have turned out to be mathematically rich objects which can be understood with the 
use of graph theoretic and linear algebraic tools [|7| . 

In section |T| we present the general construction. In section |T| we present two examples 
and discuss various open problems. 



II. UNEXTENDIBLE PRODUCT BASES AND INDECOMPOSABLE MAPS 

A 72-dimensional complex Hilbert space is denoted as 7i n . The set of all linear operators 
on a Hilbert space H„ will be denoted as B(H„). The subset of Hermitian positive semidef- 
inite operators is denoted as B(7i n ) + . We will use the conventional bra and ket notation in 
quantum mechanics, i.e. a vector ip in 7i n is written as a ket, 
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m e H n , (i) 

and the Hermitian conjugate of r/>, ip*, is denoted as a bra The complex innerproduct 
between vectors \ip) and \<p) in Ti n is denoted as 

(V#> = (2) 
The vectors |^) G are usually normalized, (?/>|V ; ) = 1- Elements of B(Ti. n ) + are denoted 



as 



p = Y^Xi\il>i){il>i\, (3) 

i 

where are the normalized eigenvectors of p and A, > are the eigenvalues. When p 
has trace equal to one, p is said to be a density matrix. The physical state of a quantum 
mechanical system is given by its density matrix. If a density matrix p has rank one, p is 
called a pure state and can be written as 

P=IM# (4) 

Let S: B(TC n ) B(H m ) be a linear map. 5 is positive when S: B(H n ) + — > B{7i m ) + . 
Let idfc be the identity map on B(TCk)- We define the map id*, Cg> S: B(Hk ® 7^ n ) — > 
B{H k ®H m ) for A; = 1,2,... by 



(id fc <g> 5) ^ (Ti <g> T^j = °i ® 5 ( r i)> 



(5) 



where <7j G B{7i k ) and Tj G B(H n ). The map 5 is ^-positive when id*, ® 5 is positive. The 
map iS is completely positive when S is /c-positive for all k = 1, 2, . . .. Following Lindblad 
||, the set of physical operations on a density matrix p G B{7i n ) + is given by the set of 
completely positive trace-preserving maps S\ B(TC n ) — > B(TC m ). Similarly as fc-positive, one 
can define a fc-copositive map. Let T:B{H n ) — > B{7i n ) be defined as matrix transposition 
in a chosen basis for TC n , i.e. 

(T(A))ij = A jh (6) 



on a matrix A e B(7i n ). The map S is fc-copositive when id^® (SoT) is positive. A positive 
linear map S: B(TC n ) — > B(H m ) is decomposable if it can be written as 

S = Si + S 2 o T, (7) 

where «Si: B(H n ) — > B(TC m ) and 5 2 : B(H n ) — > B{7i m ) are completely positive maps. It 
has been shown by Woronowicz that all positive linear maps 5: B(7i. 2 ) ~^ B(7i 2 ) and 
5: B(y-t 2 ) — > B{7i^) are decomposable. 

Definition 1 Lei p be a density matrix on a finite dimensional Hilbert space Ha ® Hb- A 
state \if)) of the form \ip A ) ® \?p B ) is a (pure) product state in Ha ® Wb. TTie density matrix 
p is entangled iff' p cannot be written as a nonnegative combination of pure product states, 
i.e. there does not exist an ensemble {pi > 0, \ipf ® ipf)} such that 

P = Eftl^X^I®l^f><v^l- (8) 

i 

When p is not entangled p is called separable. 

The problem of deciding whether a bipartite density matrix p on Ha ® is entangled 
can be quite hard. The following theorem by the Horodeckis @ formulates a necessary and 
sufficient condition for a density matrix p to be entangled: 

Theorem 1 (Horodecki) A density matrix p on Ha <8> Hb is entangled iff there exists a 
positive linear map S: Hb — > Ha such that 

(id A ® S) (p) (9) 

is not positive semidefinite. Here i&a denotes the identity map on B(Ha)- 

Remark An equivalent statement as Theorem [I] holds for positive linear maps S: Ha — * 
Hb and the positivity of S £g> ids . 

The consequences of Theorem [I] and Woronowicz' result is that a bipartite density matrix 
p on H 2 ®H 2 and H 2 ®H% is entangled iff (id^ ® <Si + 5 2 o T) (p) is not positive semidefinite. 
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As Si and S2 are completely positive maps this is equivalent to saying that (icU ® T) (p) is 
not positive semidefinite. 

The more complicated structure of the positive linear maps in higher dimensional matrix 
algebras, namely the existence of indecomposable positive maps, is reflected by the existence 
of entangled density matrices p on TCa®'Hb for which (id^ <8> T) (p) is positive semidefmite. 

The first example of such a density matrix on Ti.2 <8> 7^4 and 7^3 Cg> was found by P. 
Horodecki . In Ref . || a method was discovered to construct entangled density matrices p 



with positive semidefmite (id.A <8> T) (p) in various dimensions dim7^A > 2 and dim 7^ > 2. 
The construction was based on the notion of an unextendible product basis. Let us give the 
definition. 

Definition 2 Let 7i be a finite dimensional Hilbert space of the form TIa ®7~Lb- A partial 
product basis is a set S of mutually orthonormal pure product states spanning a proper sub- 
space Tis ofTC. An unextendible product basis is a partial product basis whose complementary 
subspace 7i§ contains no product state. 

Remark This definition can be extended to product bases in Ti = "Hi with arbitary 
m. Note we restrict ourselves to sets orthonormal sets S. 

With this notion we can construct the following density matrix: 

~ I S I 

Theorem 2 []5[] Let S be a bipartite unextendible product basis {\oij) (g) |/3i)}|Jj in Ti. = 
Ha ® T~Lb- We define a density matrix p as 

P = dim -^_ jg| f id AB - E ® IA)(A|^ , (10) 

where id^^ is the identity on Ti. The density matrix p is entangled. Furthermore, the state 
(idyi <8>[Sx+T o<S 2 ])(p) for all completely positive maps S\ and S2, is positive semidefinite. 

Proof The density matrix p is proportional to the projector on the complementary sub- 
space 7Yg . As S is unextendible contains no product states. Therefore the density matrix 
is entangled. It is not hard to see that (id^Cg>T)(p) is positive semidefinite. It has been proved 



in Ref. |TT| that when (icU ® T) (p) is positive semidefinite then (icU®T o5 2 )(p) where S2 is 



any completely positive map, is also positive semidefinite. Therefore (id/i® [«Si + T o S 2 ])(p) 
is also positive semidefinite. □ 

We are now ready to present our results relating these density matrices obtained from 
the construction in Theorem |2| to indecomposable positive linear maps. We will need the 
definition of a maximally entangled pure state: 

Definition 3 Let Ti = T~Ca <E> 7~Lb- Let be a normalized state in 7i and 

p^ = Tr B |^|, (11) 

where Tr# indicates that the trace is taken with respect to Hilbert space Tis only. The state 
ETi is maximally entangled when 

S{pA,i>) = -Tr p A ,v !og2 PA,i> = log 2 min(dimHA,dim^ i j) (12) 

The function S(pA,ip) is the von Neumann entropy of the density matrix pa,^>- 

Remarks For pure states the von Neumann entropy of p^ is always less than 
or equal to d = log 2 min(dim7-d, dim7^B). For maximally entangled states we will have 
PA,ip = diag(l/<i, . . . , 1/d, 0, . . . , 0) so that the maximum von Neumann entropy, Eq. fll2l) , is 
achieved. When dim 7-^4 = dim7is one can always make an orthonormal basis for 7i with 



maximally entangled states [12 



The following lemma bounds the innerproduct between a maximally entangled state and 
any product state. 

Lemma 1 Let 7i = H.a ®1~tB- Let G Ti be a maximally entangled state. Let d = 
min(dim?i J 4, dimHs)- For all (normalized) product states \4>a) ® \4>b), 

IWa)®|<MI 2 <^. (13) 
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Proof We write the maximally entangled state in the Schmidt polar form [|TJ] as 

m = ^j:\a t )^\b t ), (14) 

where (a,i\aj) = and (bi\bj) = 5ij. Thus we can write 

d 



2 1 



|(*|0 A )®|0 B )| ^ 



i=l 



4 (l5) 



using the Schwarz inequality and Ylf=i I (0a I a i) 1 2 < 1 an d |(0s|^)| 2 < 1- 1=1 
We will also need the following lemma: 

ISI 

Lemma 2 Let S be an unextendible product basis {\ati) ® |/3j)}i=i inH = Ha <8> "Wb- -^ei 

/(i0A),l^)) = EK^I«i)l 2 K^IA)l 2 . (i6) 
i=i 

T7ie minimum of f over all pure states \<j>A) € Ha and \4>b) £ Hb exists and is strictly 
larger than 0. 

Proof The set of all pure product states \4>a) ® \(/>b) on H is a compact set. The function 
/ is a continuous function on this set. Therefore, if there exists a set of states \4>a) ® \^>b) 
for which / is arbitrary small then there would also exist a pair \<p' A ) <S> \4>' B ) for which / = 0. 
This contradicts the fact that S is an unextendible product basis. □ 

The following two theorems contain the main result of this paper. 

I Si 

Theorem 3 Let S be an unextendible product basis <S> |A)}i=i inH = Ha ®Hb- Let 
p be the density matrix 
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|S| 



p = dimH _ |g| ( id ^ - E k)H ® IAXAI ) , (17) 



Let d = min(dim Ha, dim Hb)- Let H be a Hermitian operator given by 

|s| 

H = E N>H ® IAXAI - <W(* I, (18) 

i=l 

where is a maximally entangled state such that 



(¥|p|¥)>0, (19) 

and 

|S| 

e= min £ \{^i)?\{M^)W (20) 

where the minimum is taken over all pure states \$a) £ TYa |<Ab) £ Hb- For any 
unextendible product basis S it is possible to find a maximally entangled state such that 
Eq. (173) holds. H has the following properties: 



TrHp<0, (21) 
and for all product states \4>a) <8> \4>b) £ 

TrH|<^)(</u|® \<f> B ){<f> B \ >0. (22) 



Proof Eq. (P2|) follows from the definition of e, Eq. (|20|) , and Lemma y. Consider Eq. 



([2~T|). As the density matrix p is proportional to the projector on 7ig , one has 

TrHp = -de(tf|p|tf), (23) 

which is strictly smaller than zero by Lemma |^ and the choice of the maximally entangled 
state, Eq. (|j"9|) . When dim 7-^ = dim7^ B there exist bases of maximally entangled states 
and thus there will be a basis vector for which (\l/|p|\l/) is nonzero. In case, say, 
dim > dimTCs, the maximally entangled states form bases of subspaces TC' = 7i' A ®H.b 
with 7i' A C Ha and dim 7-^ = d\m7i B . This completes the proof. □ 

ISI 

Theorem 4 Let S be an unextendible product basis {\a>i) ® |A)}i=i inTi = Ha ® 'Hb- Let 
H 6e defined as in Theorem^ Eq. (|T^). Choose an orthonormal basis {\i)}f™ nA for Ha- 
Let S: B(TCa) ~^ B(H,b) be a linear map defined by 

S(\i)(j\) = (i\K\j)- (24) 
Then S is positive but not completely positive. S is indecomposable. 



Proof The relation between S and H, Eq. (24[), follows from the isomorphism between 
Hermitian operators on H,a ®1~Lb with the property of Eq. ( p2|) and linear positive maps, 
see p|,p^|. In particular, iff a Hermitian H operator on Ha ®Hb has the property of Eq. 
(p2|) then the linear map 1Z: B(Ha) B(Hb) defined by 

H=E(i»*>oir®^(i<>ou (25) 

is positive for any choice of the orthonormal basis {\i)}f™' HA . The map S = 1Z o T where T 
is matrix transposition in the basis {\i)}f™ nA of Eq. ([24|) is then positive as well. 

We will show how the density matrix p derived from the unextendible product basis, 
Eq. fllTD shows that S is not completely positive. At the same time we prove that the 
assumption that S is decomposable leads to a contradiction. We can rewrite Eq. fl24}) as 



H=(id yl ®5)(|^ + )(^ + |) (26) 

where \^ + ) is equal to the (unnormalized) maximally entangled state Y^=i Ha \i) ® N)- Let 
S*: B(TCb) — > B(Ha) be the Hermitian conjugate of S. We use the definition of S* 

TrS*(A*) B = TtA* S(B), (27) 

and Eq. (p6|) to derive that Eq. (|2lD can be rewritten as 

TrHp= ($+| (id A ®5*)(p)|^+) < 0, (28) 

Thus «S* cannot be completely positive and therefore S itself is not completely positive. If S 
were decomposable, then S* would be of the form S\ + T0S2 where Si and S2 are completely 
positive maps. The density matrix p is positive semidefinite under any linear map of the 
form Si + T o S2 by Theorem 0. This is in contradiction with Eq. fl28| ) and therefore S 
cannot be decomposable. □ 

We will now show how one can determine a lower bound on the value of e , Eq. (pop. 
Note that when we determine a lower bound e > e max , then all operators H, as in Eq. 
of the form 



|S| 

H = J2 k)H ® IAXAI - (29) 
i=i 

where /i G (0, e max ] correspond to positive maps, 
is 

Let {\ati) ® l/Jj)}^! be an unextendible product basis in Ha <8> T^b with ci^ = dim 7-^4 

isi ISI 

and ds = dim7i B . Let = {_ [ Qr») }4==i and = {|A)}i=i- We pick a vector and order 
the innerproducts |(aj|^)| 2 in an increasing sequence; let us call them x\ < %2 < ■ ■ ■ < x\s\. 
Then we select vectors \oti) corresponding to the smallest innerproducts in this sequence up 
to the point where the set of selected vectors \otj) spans the full ^-dimensional Hilbert space 
Ha- Let us call this set G Sa- If we would take away anyone state from S^, the remaining 
vectors would no longer span Ha- As the vectors in the set span Ha, it must be that 
X| S P| > 0. Let us label this corresponding vector as \ai max ), i.e. X\ s p\ = \{oLi max \4> a)\ 2 ■ Now 
we construct a subset of Sb in the following way; we define = | G" S^}U{\p imax )}. 
We note that the vectors in the set span the Hilbert space Hb', if not, then there would 
exist a vector \<Pb) which is orthogonal to all vectors in and a vector \4>a) which is 
orthogonal to all vectors in \ {|cKj ma:E )}, wich would in turn imply that e = 0, in other 
words the set S would be extendible. Let us pick a vector \<Pb) and denote the innerproducts 
|(A|0s)| 2 with |/3j) G as t/i < t/ 2 < . . . < 2/|Sg[- As the vectors in span Hb, we know 
that ?/| g p| > for any state \4>b)- This implies that for a particular \4>a) and \4>b) we can 
bound 

^|( ai i0A)i 2 i(A|0 fl }| 2 >^|2/|s||, (30) 

i 

the product of the two largest innerproducts of the vectors \4>a) and \4>b) with the vectors 
from and respectively. Therefore e itself, Eq. (|20|), can be bounded as 



e > min x\ai>\y\ap\, (31) 



where X\$p\ denotes the largest innerproduct between \<Pa) and a state in the set S^. We 
minimize over \4>a)— > Sa and \4>b)— > S# where the arrow denotes that a state \4>a) gives 
rise to a set as in the construction given above. A set (and similarly S^) might 
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not be uniquely denned given the vector \4>a), for example when several innerproducts of 
\4>a) with states jaij) are identical. Since the lowerbound given in Eq. ( |30D works for all 
sets S*l and which are constructed with the method given above, we could do an extra 
maximization over S% and Sg, given the states \4>a) and \4>b), but for the sake of clarity 
this maximization is omitted in Eq. (|5TD . 

We have the following proposition that can be used to bound X| S p and y\$p\ given the 
sets S^ and S^: 

Proposition 1 Let {IV^KLi ^ e a set °f n vectors in TC such that {\^i)}2=i span the Hilbert 
space Tt. Then for any vector |0) we have 

nmax|(0|^)| 2 > El(^>| 2 > X ™im (32) 
% 

where \ m in is the smallest eigenvalue of the Hermitian matrix P = J2i I'^KV'il- 
Summarizing, we get the following 

e > mm p A p B = e max . (33) 

S B' S B A\ PBl 

In order to carry out this calculation, we first find all minimal 'full rank' subsets S^ of Sa- 
Then for each of these sets S^ we compute the smallest eigenvalue of Dies^ \ a i)( a i\- Also 
for each set S^, we construct complementary sets S^ which contain all the vectors |/3j) such 
that \cti) S^ and a single state |$) such | ck<) € S^. For each set S^ there will be |S^| of 
such sets Sg. Then for each S^ we compute the smallest eigenvalue of Si<=s| Then 
we can take the minimum over all these values to obtain a bound on e. Note that this is 
now a minimization over a discrete number of values. If the set S has few symmetries and is 
defined in a high dimensional space, the procedure will be cumbersome. In small dimensions 
or for unextendible product bases which do have many symmetries, the calculation will be 
much simpler. In the next section we present two examples of positive maps based on the 
construction in Theorem f| and for one of them we will explicitly compute a lower bound on 



€. 
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III. EXAMPLES AND DISCUSSION 



As we have shown the structure of unextendible product bases carries over to indecom- 
posable positive linear maps. In this section we will list some of the results that have been 
obtained about unextendible product bases. We will take two examples of unextendible 
product bases and demonstrate the construction of Theorem and Theorem 

1. In Ref. || it was shown that there exist no unextendible product bases in 7i 2 ®H n for 
any n > 2. 

2. In Ref. |§ it was shown how to parametrize all possible unextendible product bases in 
H.3 <8> as a six-parameter family. 

3. In Ref. a family of unextendible product bases, based on quadratic residues, in 
H-n ® Tin where n is any odd number and 2n — 1 is a prime of the form 4m + 1 has 
been found. 

4. In Ref. || a family of unextendible product bases 7i n ®7i m (m > 2, n > 2) for arbitary 
m / n as well as even n = m has been conjectured. The conjecture was proved in 

<8> T~tn and Ha <8> Ha (see also fT5|j). 



5. In Ref. it was shown that when Si and S2 are unextendible product bases on 
H\ ® H} B and 7i\ ® 7i 2 B respectively, then the tensorproduct of the two sets, Si ® S 2 , 
is again an unextendible product bases on (TC\ <8> TC\) <S> (H B ® Ti.%). 

Example 1: One of the first examples of an unextendible product basis in H3 ® H3 was 
the following set of states [|| . Consider five vectors in real three-dimensional space forming 
the apex of a regular pentagonal pyramid, the height h of the pyramid being chosen such 
that nonadjacent apex vectors are orthogonal. The vectors are 

\vi) = iV(cos— ,sin— ,h), z = 0,...,4, (34) 
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with h = ~ y 1 + y/E and AT = 2/y 5 + \/5. Then the following five states in 7^ 3 <g> 7i 3 form 
an unextendible product basis: 

\Pi) = \vi) ®\v 2imod 5), i = 0,..., A. (35) 

Let p be the entangled state derived from this unextendible product basis as in Eq. ( |T0D 
Theorem 0. We choose a maximally entangled state here named \^ + ), 

|*+> = -L(|ll> + |22> + |33». (36) 

One can easily compute that 

(*+| p = -(l- 7 + V t ) > 0. (37) 
V lHl ' A\ 3(3 + ^)/ 

Let us now compute a lower bound on e, as in Eq. (|33|). Due to the high symmetry of 



this set of states, we will only need the compute the minimum eigenvalue of the Hermitian 
matrix Pi = |t>o)(fo| + + and P 2 = \vq)(v \ + |t>i)(t>i| + |v 3 )(v 3 |; all other 

subsets of three vectors, either on Bob's or Alice's side, correspond to matrices with the 
same eigenvalues as Pi or P 2 . Easy computation shows that Pi has the smallest eigenvalue 
which is equal to 



Amira — ~ • (38) 



Then as the states on Bob's side are identical, we get 

> ^nin = A + V2-V5-VT0 

e ~ 9 9 1 ' 

The map S as defined in Eq. (^4j) Theorem f|, follows directly: 

4 

5 (l*)01) = H(^l^)Kb')l t; 2fcmod5)(^2femod5| ~ 3//|z)(j|. (40) 
fc=0 

where 

4 + y/2 - V$ - VTO 



o, 
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(41) 



A positive linear map S: B(TC n ) — > B(TC m ) is unital if <S(id n ) = id m . We will demonstrate 
that S is not unital. One can write 

4 

S(id A ) = Tr A H = J2 K. mod 5) (v 2k mod 5 | - 3/i Tr A |^+) (*+ 1 , (42) 

fc=0 



which in turn is equal to 



<S(id A ) = diag ( ^^=, ^^=, \/5^ - /iid B . (43) 



The next example is based on a more general unextendible product bases that was 
presented in Ref. ||. 

Example 2: The states of S in 7i 3 ® are 

1 n— 1 

|F fc °> = ^=|0>®(|l> + £^ 2) IO), l<fc<n-3, (44) 
|F fc 1 ) = ^i=|l)®(|2) + X:^- 2) IO) ) l<A;<n-3, (45) 



1=3 

n-1 



I^ 2 ) = ^=^|2)®(|0) + E^ 2) I0), l<*<n-3, (46) 
Vn-2 l=3 

l^ 3 ) = 4(|0)-|l))®|0), (47) 



hM = ^(|l>-|2»®|l>, (48) 
l^ 5 ) = ^ ! (|2)-|0))®|2), (49) 



1 2 n-l 

I^> = -s=EEI0®Ij>, (50) 

v<j^ i=o j=o 

and we have uj = exp(2iii/(n — 2)). Here the states {\k)}^ZQ form an orthonormal basis. In 
total there are 3n — 5 states in the basis. We choose a maximally entangled state, again we 
take \^ + ), Eq. fl3"6"|). One can show that 

^+\ p \y+) = Ul--L) >0 . (51) 



5 V2 3rt 
The map 5: B{H^) — > B{7i n ) is given as 

n-3 2 6 

5(|<>01) = E E(«>(^IJ> + E(#P>(^li) - e 1001- (52) 

fc=l p=0 p=3 
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The following questions concerning the positive maps that were introduced in this paper 
are left open. 

1. Is S always non unital? We conjecture it is. As we showed, see Eq. P2|), the answer 
to this question depends on whether 

|S| 

EIAXAI ocid B , (53) 
i=i 

ISI 

where the set of states {|/3i)}|=i are one side of the unextendible product basis. The 
states will span 7i B but they will not be all orthogonal, nor all nonorthogonal. 

2. It was shown in Theorem |] that the new indecomposable positive linear maps 
S: B(H m ) — > B{7i n ) are not m-positive, as they are not completely positive. Are 
these maps S A;-positive with 1 < k < ml The answer to this question will rely on a 
better understanding of the structure of unextendible product bases. 

3. In H a single example was given of a entangled density matrix on 7i 3 ® 7i 4 which 
was positive semidefmite under ids ® T. The density matrix was based not on an 
unextendible product basis, but a 'strongly uncompletable' product basis S. It could 
be shown that the Hilbert space Hg had a product state deficit, i.e. the number of 
product states in 7Yg was less than dim7ig. It is open question how to generalize this 
example and whether these kinds of density matrices will give rise to more general 
family of indecomposable positive linear maps. 

The author would like to thank David DiVincenzo for fruitful discussions, Alan Hoffman 
for suggestions to improve the text and Peter Shor for sharing his insights in how to derive 
a lower bound on e. 
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